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Abstract 

> . 

(^T) , Using a quantumlike description for light propagation in nonhomo- 

CN ■ geneous optical fibers, quantum information processing can be imple- 

mented by optical means. Quantum-like bits (qulbits) are associated 
to light modes in the optical fiber and quantum gates to segments of 
the fiber providing an unitary transformation of the mode structure 
along a space direction. Simulation of nonlinear quantum effects is 
^ also discussed. 

£i : PACS: 03.67.Lx 

a 

1 Introduction 

> ■ 

^ ■ The quantum computer idea |], |[ Q uses the possibility to code numerical in- 

formation by vectors in Hilbert space. In the simplest case, a two-dimensional 
Hilbert space is said to code a qubit. A physical realization of a qubit might 
be a spin 1/2 particle or a two- level atom. Calculations are carried out by 
unitary time-evolution transforming an input state vector into a final state 
vector in Hilbert space. Several characteristic features distinguish the quan- 
tum computation paradigm from classical computation. First, because of 
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the superposition principle, the qubit space contains all the complex linear 
combinations, rather than just two states as in a classical bit. Second, in a 
quantum space of n qubits (of dimensionality 2 n ) there are both factorized 
and entangled states, the latter having no correspondence in the space of n 
classical bits. Third, quantum evolution operates simultaneously in the 2™ 
states, implying intrinsic exponential parallelism of quantum computation. 
To extract the result from a quantum computation one has to observe the 
system, that is, to project it in one of the exponentially many states, thus los- 
ing most of the exponential amount of information generated by the unitary 
evolution. However, it is possible to take advantage of the exponentiality 
of quantum computation using the interference mechanism characteristic of 
quantum mechanics. In short, it is the combined effect of entanglement, ex- 
ponential parallelism and interference that may allow the full potential of 
quantum computation to be realized M. 

According to the Church thesis, a classical Turing machine can simulate 
the computation of any computable function. Therefore a classical computer 
can simulate any computation of a quantum computer. The problem is how 
long the simulation will take to run. In particular it is know that some 
problems that are computed in exponential time by classical computers might 
be solved in polynomial time by quantum computers. Therefore, to take full 
advantage of the quantum computation algorithms, one requires physical 
devices obeying quantum laws. 

There are many classical systems that physically implement some of the 
features of quantum computation [|, [| ^ §]. For example electromagnetic 
waves may be linearly superimposed and interfere. Nevertheless, given the 
qualitative differences between classical and quantum mechanics, no classical 
system, where computations correspond to evolution in real time, may ever 
implement simultaneously all the features of quantum computation. Other- 
wise we would have proved the physical equivalence of classical and quantum 
mechanics. 

Notice however that quantum computation is not quantum mechanics. 
Quantum computation is a mathematical algorithm that uses all the math- 
ematical features of quantum mechanics. In particular it is irrelevant for 
the algorithm how the Hilbert space is physically implemented and whether 
the unitary evolution is taking place along a real time direction or along 
some other coordinate. It is here that quantum-like systems may play a role. 
Quantum-like systems are classical systems which obey equations formally 
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identical to the Schrodinger equation, but where the role of time is played by 
a space coordinate. Therefore, insofar as they obey equations mathematically 
identical to those of quantum mechanics, they may implement all the quan- 
tum computation operations, provided the unitary evolution is interpreted 
not as evolution in time, but as evolution along a space coordinate. Because 
of this exchange of the role of the coordinates, there is no contradiction with 
the non-equivalence of classical and quantum mechanics. 

In Sect. 2 we give an overview of quantum-like systems and then, con- 
centrating on fiber optics phenomena, we discuss how quantum information 
may be coded on the fiber and the kind of non-homogeneities and interactions 
needed to implement a set of universal quantum gates. Here we have concen- 
trated on harmonic light modes. An alternative scheme might be developed 
based on soliton propagation on the fibers. 

The fact that the unitary evolution needed for a quantum gate is obtained 
by setting up a nonuniform refraction index profile along the fiber, means 
that the unitary evolution becomes permanently coded in the hardware. This 
may be physically more convenient than to control a sequence of operations 
each time the gate is operated. 

Abrams and Lloyd have suggested that a stronger computation model 
would be obtained with non-linear terms in the quantum evolution. So far, all 
experimental evidence favors strict linearity of quantum mechanics. However, 
this is not so in quantum-like systems where non-linear effects may easily 
be introduced. Actual implementations of non-linear quantum computation 
might therefore use quantum-like systems. 

2 Quantum-like systems 

More than half a century ago, Fock and Leontovich have shown that paraxial 
beams of electromagnetic radiation in the parabolic approximation can be 
described by a Schrodinger-like equation |10|, |TIf . The role of time in this 
equation is played by the spatial (longitudinal) coordinate of the light beam, 
the role of Planck's constant is played by the light wavelength, and the role 
of potential energy by the index of refraction of the medium. Thus, the 
paraxial beam of light, a purely classical object, is a quantum-like system 
obeying equations formally identical to those of quantum mechanics. 

Given the Helmholz equation for a component of the electric field, ob- 
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tained for a fixed frequency, neglecting media dispersion and polarization 

d 2 E d 2 E 



+ — + k 2 n 2 (x,z)E = 



d 2 x dz 2 



(1) 



(we use the planar configuration, A = 2n/k is the wavelength in vacuum, and 
z the longitudinal coordinate). 

Introduce the complex function ip(x,z), which is the slowly varying am- 
plitude of the electric field in 



E(x, z) = n 1 ^ 2 (z) ip(x, z) exp 



ik / n (0 d£ 



(2) 



The ansatz (0) reduces the Helmholz equation (ffl) to a Schrodinger-like equa- 
tion 



dip(x, z) 



A d 2 ip(x,z) 



+ U (x, z) ip(x, z) 



(3) 



dz ^txtlq^z) dx 2 

U(x, z) being an effective potential related to the index of refraction of the 
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medium, n(x, z) 



U{x,z) 



n {z) 



n (z) — n (x, z) 



and Uq{z) = n(0, z) the index of refraction at the beam axis. In deriving this 
equation, second-order z— derivatives of ip and derivatives of the function 
710(2) were neglected. This is justified for slow variation of the index of 
refraction along the beam axis over distances of order of one wavelength 
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nl(z) 



dno(z) 



dz 



«1. 



The Fock-Leontovich approximation is a basis for the description of light- 

(4) 



beam propagation in optical fibers [12, 13, 14, 15] leading to 

dip(x, y, z) A 2 
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7lo(2) 



nl(z) -n 2 (x,y,z) ip(x,y,z) 



which is Schrodinger-like for a wave function depending on the transversal 
coordinates x and y. The longitudinal z coordinate plays the role of time 
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in the Schrodinger equation. The unitary z-evolution of the electromagnetic 
complex amplitude is described by the evolution operator U(z) 



U 0, z )ip(x, y, z ) = ip(x, y, z), 



associated to the Hamiltonian 



H(z) = ^ + ^ 




+ U{x,y,z). 



(5) 



with p x = —iX ,p y = —iX ^ and a potential function 



U(x,y,z) 
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n {z) 



nl(z) -n 2 (x,y,z) 



Other quantum-like systems are reviewed in [|lg, [L7j. An important ex- 
ample is sound-wave propagation in acoustic waveguides [|18| . Acoustic waves 
in the paraxial approximation are well described by a Schrodinger- like equa- 
tion. Charged-particle beams were also recently discussed as quantumlike 
systems [fijl [2(J . Light beams inside diode lasers have been treated as quan- 



tumlike systems as well, due to the waveguide structure of their active re- 
gion |21], [22 1 . As remarked in p3| the variety of quantumlike classical systems 
provides a wide range of possibilities for perfect simulation of quantum com- 
putation operations by classical means. 



3 Modes and gates in optical fibers 

In an optical fiber, the light modes, being solutions of the Schrodinger-like 
equation (^), have all the properties of quantum-mechanical wave functions 
including the entanglement phenomenon. 

Light modes in the fiber are used to code quantum- like bits (qulbits). How 
many qulbits may live in one optical fiber? In the simplest case, which is con- 
sidered here, light modes with a fixed frequency are considered. But, in the 
same optical fiber, light of different frequencies may be used. The Helmholtz 
equation holds for each frequency, with an index of refraction profile that may 
be different for different frequencies. Hence, in the same optical fiber, one 
may store many qulbits simply by exploiting the propagation of light beams 
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with different frequencies. Interaction of photons with different frequencies 
may provide useful computation effects. 

The Fock-Leontovich approximation is obtained from the Helmholtz equa- 
tion for the components of the electric or magnetic field. This equation is a 
scalar approximation which neglects the tensorial structure of the dielectric 
constant and the polarization. Also neglected are time and space dispersion, 
related to the nonlocal linear response of the medium to electromagnetic per- 
turbations. Taking into account these effects would provide an even richer 
framework to accommodate qulbits in the fiber. 

We consider now the scalar fixed frequency situation described by Eq. (f|) . 
When the index of refraction profile has the form of an inverse well, the op- 
tical fiber traps discrete modes ip nin2 (x, y, z), rii and being integer labels. 

With light modes on a fiber and the unitary z— evolution associated to 
Eq.(Q) one may perform quantum computation over continuous variables in 
a way similar to the one proposed for time evolution in Ref . |[24|| . All the 



physical interactions needed for the construction of polynomial Hamiltoni- 
ans are available by the choice of the appropriate refraction profile and by 
Kerr interactions. Also, as explained below, by restricting oneself to finite- 
dimensional subspaces of excitations, one may perform the same operations 
as in quantum computation with discrete variables. 

There are several ways to code information by light modes in a fiber which 
may be useful for quantum information processing, in particular those asso- 
ciated to different choices of basis in self-focusing potentials. The important 
self-focusing case is associated to quadratic potentials of the form 



U(x, y, z) = a(z)x 2 + b(z)y 2 + d(z)xy + e(z)x + f(x)y + l(z), 



(6) 



and in this case, an explicit solution may be obtained for the four z— dependent 
integrals of motion |T5f 






U(z,z Q ) t \U-\z,z ), f u : =U(z,z ) f \U-\z,z ) 




Defining the boson integrals of motion 

( a ( z ) \ _ 1 ( %o( z ) +ipxo(z) \ / b(z) \ 1 / y (z) + ip y0 (z) 



(8) 
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several basis may be constructed. The discrete Fock-state modes are solutions 
to the eigenvalue equation 



a)(z)a(z) | 711,712, z) = ni | nx,n 2 ,z), (9) 
b*(z)b(z) | nx,n 2 ,z) = n 2 \ nx,n 2 ,z), n 1: n 2 = 0, 1, 2, . . . (10) 

The Fock-state modes are obtained from the fundamental mode | 0, 0, z) by 

| nx,n 2 ,z) = 1 ; | 0,0, z). (11) 

Vni!n 2 ! 

A spin-like description of the Fock-state modes is related to SU (2)-subgroup 
of the Weyl-symplectic group in two dimensions. This is related to the 
Jordan-Schwinger map 

J + (z) = a\z)b{z); J_(z) = b\z)a{z); J 3 (z) = i (a\z)a{z) - b\z)b{z)) 

(12) 

Irreducible representation spaces for this subgroup are spanned by states 
{| rii,n 2 , z); ni + n 2 = N} for each fixed N. 

On the other hand, coherent modes in the optical fiber are labelled by 
two complex numbers a and j3, 

| a, (3, z) — exp aa){z) — a*a(z) exp (3b* (z) — (3*b(z) \ 0,0,2), (13) 

Using this basis, the self-focusing fiber could be considered a Gaussian chan- 
nel for numerical information. For arbitrary choices of the complex numbers 
a and (3 this is an overcomplete set. However choosing the numbers on the 
von Neumann lattice [f26[ 

Ofm.ni = ( n i + ttnmx) ; (3m 2 ,n 2 = -j= («2 + i2iTm 2 ) (14) 

and excluding two pairs (for example nix — n>x = m 2 = n 2 = 0) one obtains 
a discrete complete set of coherent modes, which might provide a basis for 
the coding of a large amount of quantum-like information in the fiber. 

We now analyze the question of what unitary transformations may be 
obtained by evolution of the light modes along the fiber. The simplest pos- 
sibility is by a change of the index of refraction profile. From Eq.M) we may 
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write a path integral representation for the evolution of the light mode along 
the fiber 

ip(x, y } z) =G (x, x , y, y , z) 4>(x , y , z) (15) 



with 

"(x,y,z) 



7rn (z) f (z)-U(t,z) 



G {x,x ,y,y ,z) = d 2 £exp{- dr 

J (x ,y ,0) { A Jo 

£ being a two-dimensional vector on the fiber sections and 



(16) 



U(X,z) = ^ n [n 2 ( Z )-n 2 (t,z)} (17) 
n \z) 

One sees from Eq.flTB]) that adjusting the index of refraction profile changes 
not only the potential but also the coefficient of the kinetic term. Let us 
consider a self-focusing quadratic potential and define, as before, a\a and 
fet, b to be creation and annihilation operators for harmonic modes along the 
x and ?/— directions respectively. For simplicity we drop the z argument in 
the operators. Then, changing the index of refraction profile gives us direct 
access to the generators 



a) a 



+ b^b\ a) + a: (a f + a) + b;(tf + b) ; (a f + a) (V + &) (18) 



By a simple reasoning using the Baker-Campbell-Hausdorff formula one con- 
cludes that by a non-uniform change of the index of refraction one may obtain 
in Eq.flTB]) all the operations of the Weyl-symplectic group in two dimensions. 
This group is not compact. Therefore the unitary representations of the full 
group are not finite-dimensional. This might be explored for information 
manipulation schemes were an unbounded number of states is manipulated. 
However for operations on a finite number of qulbits, it is the compact sub- 
groups that are important. In particular a useful subgroup is the SU (2) 
group described before in Eq.(^). Finite-dimensional irreducible spaces for 
this subgroup are 

a tn 6 tm |0>, n + m = 2k (19) 

for k = 0, |, 1, |, 2, ... with dimension 2k + 1. In this finite-dimensional spaces 
all unitary operations may be implemented on the fiber by changing the index 
of refraction profile. 
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To perform universal quantum computation it is necessary, at least, to 
have arbitrary unitary transformations on a single qulbit and a CNOT oper- 
ation on two qulbit s|F27|]. According to the discussion above the |0] , |1] qulbit 



states may be coded, for example, as follows 

|l] = a t |0>, |0]=6 f |0> (20) 

This being a k = | two-dimensional representation, all unitary transforma- 
tions may be performed in this space by the SU (2) subgroup. 

For the CNOT operation we may code the |1] —state of the control bit 
as the application of two energy quanta along the x— direction and the 
|0] —state of the control bit as the application of two energy quanta along the 
y— direction. For the target bit we use the same coding as in (EOT). Therefore 



|H] = a t3 |0 >, |01] = a^\0 >, |10] = a f 6 t2 |0 >, |00] = 6 t3 |0 > (21) 

the first label in |a/3] being the label of the target qulbit and the second the la- 
bel of the control qulbit. The subspace spanned by (|2i|) is a four-dimensional 
SU(2)— irreducible subspace. Therefore all unitary transformations may be 
implemented in this subspace and, in particular, the CNOT operation. 

We have therefore proved that, with this coding, the self- focusing fiber 
is capable of universal quantum computation. With higher order potentials 
many other possibilities would be available. For example, an interacting term 
of the form rja^ab^b appears in the Kerr Hamiltonian. This allows to perform 
phase operations on a) modes gated by the excitations. For example by 
coding the target qulbit as 

|l] = J=(l- a t)|0>, | ] = -L(l + at)|0> 

and the control qulbit as 

|l] = 6 t |0>, |0] = |0> 

the operator exp (jna^ab^bj implements the CNOT gate. 

The self-focusing fiber is a versatile medium to code quantum-like infor- 
mation and this is the reason we have emphasized the existence of several 
light mode basis. In Eq.(pOf) above, a qulbit is coded using modes in the 
x and y directions. We might as well have used a;— modes only and coded 
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the qulbit using the quadratures x and p x . Then, light propagation with the 
symmetric self-focusing Hamiltonian corresponding to the operator 



a Hadamard gate transformation. Observation of such transformation can 
be done measuring the position and direction of rays in the optical fiber for 
Gaussian wave packets. 

4 Conclusions 

1) Classical systems that are quantum-like, in the sense that their evolution 
along a space direction is described by a Schrodinger equation, possess a 
high potential for information processing including quantum computation. 
A promising system of this type consists of a light beam propagating along 
an optical fiber. It should be noticed however that similar possibilities exist 
with other systems, for example acoustic waves propagating along an acous- 
tic waveguide. Practical implementation problems to be addressed are the 
choice of the optical fiber and the definition of a coding standard for the 
qulbits. The variety of materials used and a fairly well developed optical 
fiber technology give us hope that the model Hamiltonians needed for the 
operations of quantum computing may be physically implemented in this 
medium. 

2) The fact that the unitary evolution of the quantum-like systems is 
associated to a space dimension, means that the unitary transformation is 
implemented in the hardware, rather than requiring a precise sequence of 
temporal operations. Also and this might be useful for mass production pur- 
poses, once a non-uniform refraction index profile is set up on the material, 
many different gates may be obtained simply by cutting fiber segments of 
different lengths. 

3) Finally we should also point out the easy possibility to imitate non- 
linear quantum mechanics by means of the nonlinear media response to the 
electromagnetic radiation. This might provide, as suggested in ||, even more 
powerful quantum computation algorithms. 
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